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ABSTRACT  (Unclassified) 

1 

This  report  presents  a  project  description  for  the  design  and  realization  of  a  real-time 
synthetic  aperture  radar  (SAR)  processor.  The  principles  of  SAR  processing  and  the 
state  of  the  art  are  described.  Conventional  SAR  processir  g  techniques  depend  largely  on 
FFT  transformation  techniques  and  devices.  These  offer  limited  possibilities  to  increase 
speed  and  performance  of  a  SAR  processor.  Therefore  the  design  of  new  algorithms  rather 
than  architectures  for  existing  algorithms  is  suggested.  After  that  new  architectures  will 
be  designed  dedicated  to  the  new  algorithms.  By  exploring  parallel  structures  of  the 
algorithms  and  by  using  VLSI  array  processing,  speed  and  performance  can  be  increased 
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SAMENVATTING  (Ongerubriceerd) 

Dit  rapport  bevat  een  project  beschrijving  voor  het  ontwerpen  en  realiseren  van  een  real¬ 
time  synthetic  aperture  radar  (SAR)  processor.  De  principes  van  SAR  processing  en  de 
huidige  stand  van  zaken  worden  beschreven.  Conventionale  SAR  processing  technieken  zijn 
grotendeels  gebaseerd  op  FFT-  transformatietechnieken  en  -devices.  Deze  bieden  beperkte 
mogelijkheden  om  snelheid  en  performance  van  een  SAR  processor  te  verbeteren.  Gesug- 
gereerd  wordt  om  in  de  eerste  plaats  nieuwe  a^oritmen  te  ontwerpen  in  plaats  van  archi- 
tecturen  voor  bestaande  algoritmen.  Pas  in  de  tweede  plaats  zuRen  nieuwe  architecturen 
worden  ontworpen  toegespitst  op  de  nieuwe  algorithmen.  Door  gebruik  te  maken  van  par- 
allelle  structuren  van  algoritmen  en  van  VLSI  array  processing,  kunnen  naar  verwachting 
snelheid  en  performance  aanzienl|jk  verbeterd  worden.  Het  ligt  in  de  bedoeling  dat  het 
onderzoek  zal  leiden  tot  een  promotie. 
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LIST  OF  SYMBOLS 


c:  speed  of  light  (meters/second) 

D:  physical  aperture  length  (meters) 

h;  height  of  the  aircraft  (meters) 

L:  synthetic  aperture  length  (meters) 

M:  number  of  samples  of  a  compressed  pulse 

N:  number  of  succesive  pulses  that  forms  a  synsthetic  aperture 

R:  distance  between  antenna  at  and  point  target  at  (xq,  s/b,  0)  (meters) 

R':  distance  between  antenna  at  and  point  target  at  (x,y,0),  with  (x,y,0)  nearby  {xo,yo,0) 

Ro'.  distance  between  antenna  at  to  point  target  at  (xo,  yot  0}  (meters) 

t:  time  (seconds) 

T:  interpulse  time  (seconds) 

v:  velocity  of  the  aircraft  (meters/second) 

At';  along  track  velocity  error  (meters/second) 

across  track  velocity  error  (meters/second) 

v„:  vector  containing  samples  of  the  echo  of  the  pulse  after  pulse  compresion 

x„;  x-coordinate  antenna  position  at  t„  (meters) 

a:  chirprate  (radians/second) 

0:  physical  beamwidth  (radians) 

0tynth-  synthetic  beamwidth  (radians) 
azimuth  resolution  (meters) 
hr'.  range  resolution  (meters) 

A:  wavelength  of  carrierwave  (meters) 

r:  pulse  width  (seconds) 

Tc.  pulse  width  compressed  pulse  (seconds) 

Uc-  angular  frequency  carrier  wave  (radiaas/second) 
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INTRODUCTION 


1.1  VLSI  array  processing  for  SAR 


Present-day  applications  of  digital  signal  processing  increasingly  demand  high  speed  and 
performance  processing  systems.  To  satisfy  these  demands  two  points  of  view  can  be  dis¬ 
tinguished.  The  first  one  is  to  ensure  that  a  processing  system  keeps  up  with  the  demands, 
irrespective  of  the  type  of  application.  The  second  one  is  to  decrease  the  complexity  of 
the  digital  signal  processing  algorithms,  irrespective  of  the  type  of  processing  system  that 
is  used.  It  is  obvious  that  neither  one  nor  the  other  will  give  an  optimal  performance  on 
its  own.  However,  a  combination  of  both  will  probably  do. 

A  useful  approach  in  the  design  of  digital  signal  processing  systems  for  a  dedicated  ap¬ 
plication  is  to  start  with  an  expedient  analysis  of  the  application  fields.  Hereby  one  has 
to  have  already  in  mind  some  alternative  processor  types  for  this  particular  application. 
Examples  are:  array  processors,  general  purpose  digital  signal  processors  and  application 
specific  processors.  The  question  whether  the  system  should  be  software  or  hardware  ori¬ 
ented  must  be  answered.  I.e.  the  boundaries  of  the  application  field  as  a  whole  have  to 
be  defined.  Then  an  algorithm  is  designed,  which  is  optimized  for  this  application  and  its 
field.  The  next  step  is  the  design  of  an  architecture  which  is  optimized  for  the  algorithm 
within  the  boundaries  of  the  application  field.  Then  the  sytem  is  implemented,  i.e.  the 
architecture  is  mapped  on  some  hard-  or  software  components. 

An  application  digital  signal  processing  is  synthetic  aperture  radar  (SAR)  processing. 
A  SAR  is  an  imaging  radar  on  a  moving  platform,  in  general  an  aircraft  or  satellite. 
SAR  processing  is  based  on  synthetic  geno'ating  of  an  effective  long  antenna  by  signal 
processing  rather  than  using  a  long  physical  antenna.  The  amount  of  SAR  data  which 
has  to  be  processed  to  generate  an  image  is  so  high  that  high  speed,  high  performance 
and  low  cost  procesnng  is  hardly  possible.  Due  to  the  progress  in  (VLSI)  digital  rignal 
procesnng  techniqms,  real-time  high  performance  and  low  cost  SAR  proMsmtg  systjems 
ornie  within  reach.  Because  of  thte  rMscm  the  Physics  and  E3ectronics  Laboratory,  as  a 
leading  institute  mi  radar  teehaology,  decided  to  (ksign  a  real-time  on-board  aircraft  SAR 


{tfocessor. 


Within  the  application  field  (SAR  processing)  two  mayor  requirements  can  be  distin¬ 
guished  directly.  The  processor  must  be  real-time  and  it  must  be  an  on-board  system. 
The  first  requirement  demands  high  processing  speed  and  the  second  one  demands  restric¬ 
tions  on  the  size  of  the  processor  and  its  power  consumption.  A  VLSI  type  array  processor 
is,  probably,  the  best  implementation  to  meet  these  requirements.  The  quality  of  the  im¬ 
age  must  be  at  least  such  that  the  flight  can  be  judged  as  successful  or  not.  VLSI  array 
processor  design  for  real-time  SAR  processing  is  a  relatively  new  subject.  Therefore  the 
first  step  of  the  design  trajectory  is  the  design  of  algorithms  (as  stated  above).  VLSI  atrray 
processing  requires  algorithms  that  have  low  computational  complexity  and  a  high  degree 
of  parallelism.  The  second  step  is  to  map  the  resulting  algorithms  onto  architectures. 
Requirements  for  the  architecture  are  a  high  degree  of  regularity  and  parallelism. 


1 .2  SAR  processing  for  VLSI  array  processors 


Traditionally  SAR  processing  algorithms  can  be  divided  in  two  mayor  parts:  range  com¬ 
pression  and  azimuth  compression.  Range  compression  is  a  signal  processing  approach 
that  increases  the  resolution  in  range  direction  (across  the  movement  trajectory).  The 
resolution  in  range  depends  mainly  on  the  pulse  width.  An  ideal  pulse  would  be  a  narrow 
pulse  with  high  peak  power,  such  that  a  good  signal- to- noise  {S/N)  ratio  is  achieved. 
In  practice  such  a  system  is  not  feasible.  An  alternative  is  to  transmit  a  relatively  long 
pulse  with  low  peak  power,  but  with  high  average  power.  After  receiving  the  echo  it  is 
compressed  to  a  narrow  pulse  with  high  peak  power.  Azimuth  compression  is  a  signal 
processing  approach  to  increase  the  resolution  in  azimuth  direction  (along  the  movement 
trajectory).  A  single  pulse  radar  would  have  a  resolution  in  azimuth  that  depends  on  the 
beamwidth  and  the  distance  between  antenna  and  ground.  By  processing  successive  echoes 
the  azimuth  resolution  is  increased  significantly,  and  becomes  also  independent  of  the  dis¬ 
tance  between  antenna  and  ground.  Range  and  azimuth  processing  are  mathematically 
equivalent  techniques. 

Beside  range  and  azimuth  compression  also  phase  error  compensation  is  of  great  impor¬ 
tance  in  SAR  processing.  Phase  errors  are  caused  by,  for  example,  uncontrolled  motions 
of  the  moving  platform,  atmospherical  losses,  etc..  These  dfects  cause,  among  others, 
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blurring  and  defocussing  of  the  image,  which  will  result  in  a  lower  resolution.  Phase  er¬ 
rors  can  be  estimated  and  compensated  by  several  image  processing  and  signal  processing 
techniques.  In  general  this  will  be  done  after  range  and  azimuth  comprrasion.  In  this 
report  a  few  examples  will  be  given,  but  details  and  a  broad  outline  of  possible  techniques 
can  be  found  in  the  references. 

Traditional  algorithms  for  range  and  azimuth  compression  and  phase  error  compensation 
have  high  complexity  and  use  large  storage  facilities.  This  implies  much  processing  time, 
large  processing  and  storage  equipment  and  much  overhead  for  memory  management, 
processing  control,  etc..  Untill  now  SAR  processing  is  performed  off-line  and  off-board, 
and  the  main  requirement  is  to  achieve  high  resolution.  Therefore  little  effort  is  given  to 
increase  the  performance,  for  example  by  designing  new  algorithms.  An  exception  is  made 
for  phase  error  estimation  and  compensation  algorithms,  but  the  reason  for  this  is  not  a 
desired  increase  of  perormance  but  an  increase  of  resolution. 

As  stated  before  real-time  high  performance  SAR  processing  comes  within  reach  because  of 
the  progress  in  V'LSI  array  processing  techniques.  Due  to  the  high  complexity,  traditional 
SAR  processing  algorithms,  and  thus  architectures,  are  not  suitable  for  implementation 
on  VLSI  array  processors.  Hence  new  algorithms  and  architectures  dedicated  to  VLSI 
array  processing  must  be  designed,  as  an  alternative  for  traditional  algorithms  for  range 
and  azimuth  compression  and  phase  error  compensation. 


t 

i 


1.3  Outline 


This  report  describes  the  design  project  of  a  real-time  on-board  SAR  processor.  New  al¬ 
gorithms  and  architectures  will  be  designed  to  perform  SAR  processing.  Th«e  algorithms 
and  architectures  must  be  dedicated  to  VLSI  array  processing  systems.  The  intended 
r«ult  of  the  project  will  be  the  definition,  design  and  realization  of  a  real-time  on-board 
SAR  processor.  The  necessary  research  for  this  project  wQl  be  carried  out  in  cooperation 
with  the  Ddft  University  of  Tbehnedogy,  faculty  of  Electrical  Ei^neering,  section  Network 
Theory.  It  is  the  intention  that  this  resotreh  will  also  lead  to  a  PhD  graduation. 
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The  orgzinization  of  this  report  is  as  foUows.  Chapter  2  gives  the  background  theory  of 
SAR  processing.  The  conventional  approach  of  SAR  processing  is  shortly  outlined,  and 
a  discussion  of  the  state  of  the  art  is  given.  A  short  introduction  of  the  Physics  and 
Electronics  Laboratory  and  the  position  of  this  project  in  it  is  given.  Chapter  3  defines 
the  project.  Objectives,  intended  results  and  the  requirements  are  outlined  and  the  sub- 
projects  are  defined.  The  project  stages  and  its  activities  are  defined  and  described  in 
chapter  4.  Chapter  5  gives  the  time  planning,  information  flows  and  organization  and 
defines  quality  norms.  Finally,  chapter  6  describes  the  planning  and  activities  for  the 
research,  which  has  to  lead  to  a  PhD  graduation. 
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2  BACKGROUND 

2.1  Synthetic  aperture  radar  theory 


Synthetic  aperture  radar.  A  synthetic  aperture  radar  (SAR)  is  a  side-looking  imaging 
radar  on  a  moving  platform  [1,  15,  18,  20,  26,  29].  Usually  the  platform  is  a  satellite 
(spaceborne  SAR)  or  aircraft  (airborne  SAR).  SAR  is  based  on  synthetic  generating  of 
an  effective  long  antenna  by  signal  processing  rather  than  using  a  long  physical  antenna. 
The  actual  physical  antenna  length  is  in  most  cases  relatively  small.  With  the  synthetic 
antenna  a  narrow  beam  in  azimuth  can  be  achieved,  and  thus  a  higher  resolution,  relative 
to  the  actual  physical  antenna.  Let  the  physical  antenna  length  be  D  and  the  wavelength 
of  the  transmitted  signal  A  than  the  physical  half-power  beamwidth  in  radians  is  given  by: 

(2.1) 

If  L  is  the  length  of  the  synthetic  array,  with  L  ">  D,  then  the  beam- width  of  the  synthetic 
antenna  is  given  by  [15,  26,  29] 


/?.ynt  =  ^  (2.2) 

Let  the  distance  between  the  antenna  and  a  point  on  the  ground  be  R.  Then  the  maxi¬ 
mum  synthetic  aperture  length  is  limited  by  the  length  of  the  illumination  pattern  of  the 
physical  beam  on  the  ground  [15,  26] 


L< 


RX 

D 


(2.3) 


This  gives  for  the  maximum  resolution  in  azimuth  (parallel  to  the  flight  trajectory),  6a, 
the  expression 
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Figure  2.1:  Geometry  of  the  SAR  system. 


Sa  —  (2.4) 

while  for  the  real  aperture  the  maximum  resolution  would  be  0R  =  XRID.  For  example 
in  airborne  systems,  with  (in  meters)  R  of  order  10®,  D  of  order  10"^  and  A  of  order  10~®, 
we  have  a  maximal  azimuth  resolution  of  order  10®  for  a  real  aperture  and  a  resolution 
of  order  10“'  for  a  sy.ithetic  aperture.  It  should  also  be  noticed  that  the  resolution  of 
a  synthetic  aperture  is  independent  of  R  and  A.  In  the  rest  of  the  document  it  will  be 
assumed  that,  if  not  stated  otherwise,  a  SAR  system  is  airborne. 

Let  the  transmitted  signal  be  a  sequence  of  pulses.  Let  p(t)  be  a  single  pulse,  defined 
in  the  interval  (-r/2,r/2)  and  zero  otherwise,  with  r  the  pulse  width.  Assume  that  at 
t  =  in,  n  =  •  •  • ,  -1, 0, 1, •  •  •,  a  single  pulse  is  transmitted.  Let  the  interpulse  time  T  be 
relatively  large:  T  r,  and  let  the  carrier  frequency  be  Uc  then  the  transmitted  signal  is 
given  as 


fl 


(2.5) 


The  geometry  of  the  SAR  system  is  shown  in  figure  2.1.  Let  a  single  unit  point  scatterer 
have  coordinates  (xq,  j/q,  0),  and  let  the  height  of  the  flight  trajectory  be  h.  Let  the  dis- 
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tanci  between  the  antenna  at  t  =  <„  and  the  point  scatterer  be  defined  as  R  and  define 
Rq  —  yjvo'^  Notice  that  R  depends  on  t„.  Furthermore  the  x-coordinates  of  the 
successive  positions  of  the  aircraft  at  t  =  are  defined  as  x  =  Xn-  Assume  that  xo  =  0 
and  let  v  be  the  velocity  of  the  aircraft,  then  x„  =  nvT.  Let  c  be  the  speed  of  light  then 
the  echo  signal  of  s{t)  can  be  described  as 


s{t - )  =  e-’  ‘)pit- in - ) 

c  ~  c 


where  2R/c  is  the  round  trip  time  delay  of  the  reflected  pulse.  Usually  the  term  is 
omitted,  since  it  contains  no  information  of  the  scatterer.  Hence 


n 


2H  2R. 

C  p(t  -U - ) 

c 


Pulse  compression.  Pulse  compression  is  a  technique  to  increase  the  resolution  in  range 
(paraUel  to  the  j/-axis).  The  resolution  in  range  is  determined  by  the  pulse  width  of  the 
transmitted  pulse.  The  resolution  in  range  without  pulse  compression  is  given  as  cr/2 
[15].  This  means  that  high  resolution  in  range  requires  a  very  narrow  pulse.  At  the  other 
hand  high  pulse  power  is  required  to  achieve  a  good  SjN  ratio.  An  ideEd  pulse  form  would 
be  an  appro.\imation  of  a  Dirac  pulse.  Due  to  the  physical  constraints  of  the  transmitter 
this  is  not  feasible.  A  solution  to  this  problem  is  to  transmit  a  linear  frequency  modulated 
pulse,  a  so-called  "chirp”,  defined  as 

p(t)  =  (2.8) 

where  a  is  the  chirprate.  Figure  2.2a  shows  a  chirp.  A  sequence  of  relative  long  pulses 
is  transmitted.  After  receiving  the  echo  it  is  subjected  to  a  matched  filter.  The  matched 
filter  compresses  the  received  pulses  to  very  narrow  pulses  with  the  same  power  as  the 
original  echo  pulse.  Let  Sc(t)  be  the  echo  signal,  subjected  to  the  matched  filter,  then 
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Figure  2.2:  (a)  Chirp  and  (b)  compressed  pulse. 


Sc(t) 


\  f  ,  2R.  2R\^ 

-  J  4(1  - (t  -  —)il 


-JWtf  ^ -  f  ^ 


r  J-r/2 


(2.9) 


where  ‘  means  the  conjugate  and  R'  is  the  distance  between  the  antenna  and  a  point 
(x,  j/,0)  nearby  (0,yoiO),  see  figure  2.1.  Expanding  the  integral  in  2.9  gives 


1 

r  J-tI2 


r  J-r/2 

=  e>“^(R*-R'*)sinc(ar-(i?  -  R'))  (2.10) 

c 


The  last  expression  is  a  compressed  pulse  with  the  same  pulse  power  as  before  the  filtering, 
but  with  much  higher  S/N  ratio  and  much  smaller  pulse  width,  see  figure  2.2b.  Then  2.9 
becomes’ 


‘In  thk  expwHioa  it  is  assunted  that  over  the  ctmiidete  ^rertnre  the  phase  change  of  2.10  is  small 
according  to  the  phase  change  w,  212^,  which  is  only  tme  for  short  aperture  lengths.  In  real  SAR 
systems  however  this  may  came  a  phM  enw  (range  migration),  but  in  this  system  the  error  is  assumed 
negketaUe. 
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a^(<)  =  (2.11) 

^  n 

The  pulse  width  of  the  compressed  pulse  is  given  by  =  2.8/or  [26].  The  range  reso¬ 
lution  Sr  is  then  given  by 


Sr  = 


CTc 

2 


1.4c 

or 


(2.12) 


Azimuth  compression.  Azimuth  compression  is  the  signal  processing  technique  that 
obtains  a  high  resolution  in  azimuth  from  a  sequence  of  echo  pulses.  From  figure  2.1,  R' 
and  R  are  approximated  as 


R'  =  ^Rl  +  ix-  x„y  ^Ro+^- 
R  =  /rI  +  xI^Ro+^ 


(2.13) 

(2.14) 


These  approximations  are  valid  if  i?  >  i  —  z„  and  R*  >  which  is  true  in  most  cases. 
Let  the  synthetic  aperture  length  be  L  then  the  number  of  pulses  transmitted  during  the 
building  of  the  aperture  is  given  as  TV  -h  1  =  L/vT.  Then  expanding  the  summation  in  2.11 
becomes 


N/2 

E 


na-N/i 


ns-JV/2 

•i«(¥Sg) 


(2.15) 
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where  x„  =  nvT  and  2uch  —  45r/A  are  used.  Figure  2.3a  shows  the  magnitude  term 
of  2.15  as  function  of  x.  The  azimuth  resolution  is  derived  from  the  half-power  bandwidth 
of  2.15  [26] 


^  1.4Ai?o 

The  maximum  resolution  is  obtained  by  substitution  the  maximum  aperture  length  L  = 
RqX! D  which  gives  6a  ^  Df2.  Notice  that  this  value  was  expected  from  2.4.  The  echo 
signal  after  pulse  and  azimuth  compression  is  thus 


'  sin(4ff??) 


sinc(ar-(i?  -  R'))  (2.17) 

c 


Figure  2.3b  shows  the  magnitude  term  of  the  complete  filtered  echo  signal  2.17  as  function 
of  X  and  R'. 


Phase  errors.  In  the  previous  it  was  assumed  that  the  flight  trajectory  was  a  straight  line, 
that  the  velocity  was  constant,  that  no  atmospherical  losses  did  occur,  etc..  In  practical 
systems  this  will  not  be  the  case,  moreover,  it  results  in  phase  errors  in  the  echo  signal. 
These  errors  cause  unwanted  effects  on  the  image,  like  defocussing.  In  [18]  phase  errors 
are  subdivided  in  two  groups;  predictable  and  unpredictable  phase  errors. 


Examples  of  predictable  phase  errors  are:  errors  caused  by  geocoding  and  range  migra¬ 
tion.  Geocoding  is  the  mapping  the  SAR  coordinate  system  (the  (r,x)  coordinate  system) 
onto  an  image  coordinate  system  (usually  the  (x,y)  coordinate  system).  An  example  of 
the  error  that  can  occur  is  the  following.  Let  R  be  the  distance  between  the  antenna 
and  a  target  on  the  ground  with  coordinates  (x,  y).  Let  a  nearby  target  with  coordinates 
(z,y  +  Ay)  have  distance  R  •¥  Sf.  The  image  plane  is  the  (x,y)  plane,  thus  Ay  is  the 
actual  resolution.  Figure  2.4a  shows  that  Ay  depends  on  the  angle  tfy,  while  6r  is  inde¬ 
pendent  of  This  dependency  results  in  a  predictable  phase  error.  More  examples  of 
errors  caused  by  ^ocoding  can  be  found  in  [18].  Range  migration  is  caused  by  a  pulse 
compression  assumption.  The  distances  between  successive  antenna  pomtkms  during  the 


Figure  2.3:  (a)  Azimuth  compression  and  (b)  complete  echo  signal. 

aperture  building  and  a  single  point  target  are  assumed  constant.  This  is  illustrated  in 
figure  2.4b.  These  distances  differs  maximally  +  K^-R.  For  long  aperture  lengths 

range  migration  causes  relevant  phase  errors.  Since  it  is  completely  predictable  however, 
it  will  be  neglected  in  this  section. 

Unpredictable  phase  errors  are  caused  by:  unknown  ground  height  variation,  atmospheric 
propagation  fluctuations  and  unknown  flight  trajectory.  From  these  three  the  last  one 
causes  the  most  important  (unwanted)  phase  error.  An  error  in  the  flight  trajectory  will 
hardly  affect  the  phase  during  the  round-trip  time  of  a  single  echo  pulse.  It  will  only 
affect  the  phase  during  several  pulses.  Figure  2.5  shows  that  the  displacement  in  the  flight 
trajectory  can  be  expressed  in  range  difference  as 


R-Ro^yit)+^  (2.18) 

Then  the  phase  shift  can  be  written  as 

v»(f)  =  x(iK0  +  |^)  (2.19) 


If  the  displacements  are  expanded  about  1  =  0  then 
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Figure  2.4:  (a)  Geocoding  phase  error  and  (b)  range  migration. 

i(l)  =  (t,  + +  •••  (2.20) 

5(1)  =  +  +  +  (2,21) 

J  j2  j  j2 

where  Ar,  are  along-track  velocity  errors  and  v,,  ,•••  are 

across-track  velocity  errors.  Hence,  with  the  assumption  that  the  terms  of  fourth-order 
are  neglectable,  the  phase  shift  becomes 

4ir 

¥>(0  =  y^  +  V>e(0  (2.22) 

where 

v»e(0  =  -h  (2.23) 

Notice  that  if  Av,  v,  Vr,  0  then  ip,{t)  =  0  and  ^t)  equals  ^(R  -  R')  in 

eqnatMn  2.9,  if  R  and  A'  are  approximated  by  2.13  and  2.14  and  ««  s  0  and  a  s  «<. 
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Figure  2.5:  Distortion  in  flight  trajectory. 

j  j 

If  the  parameters  Av,  ^v,  Vr,  ^Vr  are  known,  for  example  by  motion  measure¬ 

ments,  it  is  possible  to  compensate  the  phase  error  V5e(0-  This  technique  is  called  motion 
compensation  and  is  discussed  in  [18].  Estimation  of  the  phase  error  from  the  blurred 
image  is  called  autofocussing.  A  great  number  of  estimation  techniques  is  known  from  the 
literature  [2,  3,  4,  5,  8,  10,  12,  18,  19,  21,  22,  25,  30).  It  goes  beyond  the  scope  of  this 
document  to  describe  them  all,  but  two  main  trends  in  autofocussing  will  be  mentioned. 
The  first  one  is  the  mapdrift  method  [5,  12].  The  aperture  is  sub-divided  in  several  sub¬ 
apertures  and  a  multiple  look  at  the  same  im^e  is  performed.  Due  to  the  phase  errors 
each  map  will  be  shifted  relatively  from  the  others.  This  phase  error  is  determined  by 
cross  correlating  the  maps.  The  order  of  the  phase  error  that  can  be  compensated  is 
determined  by  the  number  of  sub-apertures.  A  disadvantage  of  this  method  is  that  the 
azimuth  resolution  decreases  inversely  proportional  to  the  number  of  sub-apertures.  The 
second  trend  is  based  on  contrast  optimization  [5,  8,  12,  21].  Methods  that  use  contrast 
optimization  are  based  on  determination  of  some  strongest  targets.  A  disadvantage  is  that 
high-clutter  or  less-contrast  images  decrease  the  performance  of  this  method. 
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Range  reference  function  Azimuth  reference  function 


Figure  2.6:  SAR  processing  in  FFT  domain. 


2.2  State  of  the  art 


Real-time  processing.  Real-time^  SAR  processors  were  first  developed  in  the  mid  1970’s 
[16],  as  a  result  of  the  advancement  in  digital  signal  processing.  The  techniques  that  were 
used  until  now  remianed  more  or  less  unchanged.  In  most  cases  pulse  compression  and 
azimuth  compression  are  performed  in  the  frequency  domain,  as  is  shown  in  figure  2.6. 
Proposed  architectures  are  based  on  FFT  processing  [6,  9,  11,  13,  23,  24,  27,  31).  Due  to 
the  progress  in  the  field  of  VLSI  for  FFT  processing  and  array  processing  the  performance 
of  the  architectures  increased.  A  few  proposed  architectures  use  other  transformation 
techniques,  like  polynomial  and  Fermat  number  transforms  [7,  28).  Implementation  of  the 
architectures  have  resulted  in  parallel  processors  and  computer  systems. 

The  computational  complexity  of  architectures  based  on  FFT  processing  is  still  too  high 
for  direct  VLSI  implementation.  To  illustrate  this  an  example  is  pven.  Successively  single 
pulse  echoes  are  sampled  and  pulse  compressed.  Let  n,  n  =  •  ■  • ,  — 1, 0, 1,  •  •  ■  denote  the 
transmitted  and  received  pulse.  Let  the  M  samples  of  echo  n  (after  pulse  compression) 
be  stored  in  a  vector  v„.  Let  the  image  line  I  refer  to  the  pixel  values  of  the  points 
m  Azimuth  compression  to  obtain  image  line  I  requires  N  +  I 

pulse  echoes  V|_„,  n  =  -iV/2,  ”',1^/2.  The  next  image  line  /  -J- 1  requires  the  N  pulse 
echoes  vj_„,  n  =  -JV/2+  l,”‘,Ar/2  phu  the  echo  of  the  next  pulse  vn/w-  Hence  a 
continuous  storage  of  fif  +  1  vectm  is  required  to  perform  azimuth  compresoon,  which 

*R«al-tiaM  pinriwiig  can  be  pecfonued  ou-bosid  or  off-bosid.  Si*oe  eflidait  opentkmsl  icsl-tiaw 
oa-bosid  SAR  pinrssrfsg  sytteBn  bsve  aot  been  praeated  yet  ia  the  fitentaie,  both  aspects  of  laal-Uiae 
pwwsssiag  sbaB  becBscasssd. 
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corresponds  to  storage  of  iV  +  1  x  Af  samples  in  the  corner  turning  memory,  with  N  and 
M  both  of  order  10^.  Note  that  within  the  interpulse  time  (of  order  10”^)  this  data  has 
to  be  shifted,  loaded  from  memory  to  processor  and  processed.  Since  processing  is  usually 
done  in  frequency  domain  this  implies  an  order  of  M{N  +  l)log(./V  +  1)  computations 
per  image  line  necessary  to  perform  the  FFT’s  (note  that  pulse  compression  requires  an 
additional  MlogM  computations). 

Real-time  systems.  Three  representative  examples  of  real-time  SAR  systems  (ordered 
chronologicalUy  on  date  of  publication)  are: 


1.  (1982)  A  real-time  spaceborne  SAR  system  for  off-board  SAR  procesing  is  presented 
in  [31].  It  is  used  to  process  raw  SAR  data  from  the  SEAS  AT  satellite.  The  system 
configuration  consists  of  a  minicomputer  that  acts  as  a  host  for  three  parallel  array 
processors,  two  disks  for  data  storage  and  one  disk  for  software  storage.  The  system 
is  software  oriented.  The  processing  time  for  a  single-look  image®  (1.2  x  10®  pixels) 
is  90  seconds  and  the  throughput  rate  is  6700  pixels  per  second. 

2.  (1987)  A  modular  reconfigurable  real-time  spaceborne  SAR  system  is  proposed  in 
[9].  It  is  tested  on  raw  SAR  data  from  the  SEAS  AT  satellite.  The  system  consists 
of  three  devices:  a  corner  turning  memory  (DRAM’s),  an  FFT  processor  (2/i  CMOS 
VLSI  technology)  and  a  systolic  processor.  The  processing  is  broken  up  in  five 
steps:  1.  FFT  of  echo  data,  2.  pulse  compression,  IFFT  and  storage  in  corner 
turning  memory,  3.  FFT  of  azimuth  data,  4.  phase  correction,  azimuth  compression 
and  IFFT,  5.  incoherently  summing  of  multiple  look-data.  Each  step  requires  a 
recirculation  through  the  SAR  processor:  from  memory  through  systolic  processor 
and  FFT  processor  and  back  to  memory.  The  processing  of  16  x  10®  samples  of 
raw  data  into  one  four-look  image  takes  2  seconds.  The  proposed  SAR  processor 
architecture  is  an  ofF-board  processor. 

3.  (1989)  A  real-time  on-board  SAR  system  is  proposed  in  [23].  A  mesh  connected 
SIMD  architecture  is  nsed  to  meet  the  real-time  requirement.  The  configuration 

*A  lo<A  M  u  imafc  that  icsalts  froei  the  astanath  pforwihig  of  the  data  <d  a  (aab-)i9ertue.  A  lia^ 
look  iaage  is  the  teeah  of  aaimth  prnrMriag  of  the  complete  apertam  data.  A  malt^ie-look  image  ie 
the  femh  of  the  a^ath  profieiiag  of  the  data  of  arah^  oah-apertaiet.  The  ieealtiii(  hnagee  of  the 
•ab-apettaieo  am  raauned  iacoheieally  to  obtain  the  oomflete  inu^. 
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consists  of  a  host,  a  processor  array,  a  data  memory,  a  program  memory  and  a 
processor  array  control  unit.  The  processor  array  contains  64  X  64  mesh  connected 
custom  processor  elements  (1,2/i  CMOS  VLSI  technology).  The  processing  of  16x10® 
samples  of  raw  data  into  one  single-look  image  takes  1  second.  The  system  follows  a 
maximal  pulse  repetition  frequency  of  lOOOHz.  The  implementation  has  resulted  in 
a  prototype  on-board  SAR  processor.  The  performance  meets  the  real-time  require¬ 
ment  theoretically,  but  power  consumption  and  small-size  hardware  implementations 
have  not  been  met  yet.  An  attempt  will  be  made  to  implement  the  processor  array 
on  a  stack  of  Interconnected  silicon  wafers. 

Discussion.  From  the  examples  follows  that  the  trends  in  real-time  SAR  processing 
are:  increasing  the  processing  speed  (from  ’’almost”  real-time  to  real-time)  and  decreasing 
the  size  of  the  system  (from  software  to  hardware  oriented  systems).  The  only  serious 
attempt  to  realize  an  on-board  real-time  SAR  processor  that  is  found  in  the  literature  is 
outlined  in  the  last  example,  although  it  has  not  yet  resulted  in  an  operational  system. 
Processing  in  FFT  domain  however  might  not  be  the  most  efficient  solution  to  reduce 
the  computational  complexity  and  the  use  of  memory.  It  is  remarkable  that  hardly  any 
attempt  is  made  to  find  alternative  approaches  for  processing  in  FFT  domain.  Time 
domain  solutions  may  be  more  sophisticated  with  the  present-day  techniques  for  parallel 
algorithm  design.  Therefore  a  new  direction  in  real-time  SAR  processing  can  be  the 
development  of  alternative  algorithms  for  SAR  processing.  More  significant  performance 
improvement  would  be  obtained  by  the  design  of  algorithms  which  are  optimized  in  terms 
of  parallelization  and  pipelinabillity.  The  basis  for  this  project  is  therefore  to  design  new 
algorithms  for  SAR  processing. 


2.3  FYamework  within  TNO-FEL 


Physics  and  Electronics  Laboratory.  The  Physics  and  Electronics  Laboratory  (FEL) 
is  one  of  the  three  institntes  of  the  TNO  Division  for  Defense  Research  in  the  Nether¬ 
lands.  Most  of  its  research  and  devdopment  is  sponsored  by  the  Ministry  of  Defense,  bnt 
TNO-FEL  also  teams  np  with  industry  in  national  and  international  defence  related  and 
dvil  projects.  TNO-FEL  is  orgauied  into  four  research  diddons  and  one  dividon  for 
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technological  developments.  These  are: 

1.  Operational  Research 

2.  System  Development  and  Information  Technology 

3.  Radar  and  Communications 

4.  Physics  and  Acoustics 

5.  Technological  Developments 

PHARUS.  TNO-FEL,  the  National  Aerospace  Laboratory  (NLR)  and  the  Delft  Univer¬ 
sity  of  Technology  (faculty  of  Electrical  Engineering,  section  Microwave  Technology  and 
Telecommunication)  started  several  years  ago  a  plan  for  the  development  of  a  C-band  air¬ 
borne  SAR,  the  so-called  PHARUS  (PHased  ARray  Universal  Sar)  [14].  The  final  goal  of 
the  project  is  a  full  polarimetric  C-band  (A  =  0.056m)  airborne  SAR.  The  system  wiU  be 
ready  in  1994.  It  will  possess  an  active  phased  array  patch  antenna,  which  enables  flexibil¬ 
ity  in  imaging  modes  in  addition  to  the  ability  to  choose  the  resolution.  Furthermore,  the 
antenna  pointing  on-board  the  moving  aircraft  is  maintained  through  the  phased  array, 
rather  than  by  mechanical  stabilization.  The  system  that  is  being  developed  will  be  used 
to  investigate  military  and  civil  applications,  mainly  in  the  field  of  remote  sensing  (earth 
observation),  PHARUS  is  a  system  that  requires  off-line  SAR  processing.  Its  power  is  that 
it  can  transmit  and  receive  two  orthogonal  polarizations.  Since  the  brightness  of  reflectors 
in  a  radar  image  depends  on  polarization,  the  contrast  can  be  optimized  and  reflectors 
may  be  classified  or  even  identified. 

Before  the  PHARUS  system  will  be  constructed,  experience  in  the  field  of  SAR  systems 
was  gained  by  the  development  of  a  limited  prototype  system,  called  PHARS,  meant  as 
a  testbed.  The  PHARS  is  not  polarimetric.  This  system  was  completed  recently.  The 
PHARS  records  a  swath  width  (image  strip)  of  7  km  at  a  maximum  range  of  13  km.  The 
azimuth  resolution  is  1.2  meter  (single  look).  The  range  resolution  is  4.8  meter.  The 
system  is  flown  on  a  Swearingen  Metro,  NLR’s  laboratory  aircraft,  at  altitudes  up  to  6  km 
and  speeds  up  to  120  m/s.  The  incident  angle  can  be  varied  between  20  and  60  degrees. 
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To  keep  in  touch  with  on-line  SAR  processing  TNO-FEL  started  with  the  development 
of  an  on-line  real-time  SAR  processor.  The  purpose  of  this  processor  is  to  display  on-line 
an  image  of  the  illuminated  strip.  The  SAR  processor  will  be  used  for  general  purpose 
stripmap  SAR  processing.  The  PHARS  is  used  as  a  testbed  for  the  SAR  processor.  In  the 
future  the  experience  gathered  with  the  SAR  processor  will  be  used  to  develop  a  complete 
real-time  on-board  polarimetric  SAR  system. 


4 


Page 

25 


3  PROJECT  DESCRIPTION 

3.1  Objectives 

Main  objective:  The  definition,  design  and  realization  of  an  on-board  real-time  synthetic 
aperture  processor  within  the  present-day  or  near-future  state  of  the  technique,  especially 
digital  signal  processing  and  VLSI  array  processing. 

Sub-objectives:  The  main  objective  can  be  refined  in  the  foUowing  sub-objectives: 

•  To  decrease  the  costs  of  SAR  data  acquisition  such  that  a  real-time  on-board  SAR 
processor  can  be  realized. 

•  To  combine  expertise  on  the  fields  of  signal  processing,  signal  processing  applications 
(radar)  and  VLSI  array  processing  such  that  the  leading  positions  in  these  fields  will 
be  maintained. 

•  To  obtain  a  PhD  graduation. 


3.2  Results 


The  project  has  the  following  intended  results: 


•  A  real-time  SAR  processor  definition  is  given.  The  definition  should  be  as  broad  as 
possible.  Restriction  due  to  hardware  limitations  are  omitted. 

•  A  real-time  on-board  SAR  processor  is  designed  and  realized.  This  processor  should 
be  an  implementation  derived  from  the  definition.  The  processor  displays  a  "quick- 
look”  SAR  image.  Quick-k>ok  means  that  quality  and  resdution  are  of  less  impor¬ 
tance.  A  global  view  of  the  illuminated  area  is  displayed.  The  equipment  can  be 
checked  on-line  whether  it  functions  wdl  or  not.  The  SAR  data  can  be  judged  on 
its  usefulness  by  using  the  information  of  the  quick-look  SAR  image. 
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•  SAR  processing  as  an  appUcation  of  digital  signal  processing  VLSI  array  pro¬ 
cessing  is  better  understood.  The  project  increases  insight  in: 

1 .  problems  that  occur  in  applications  of  digital  signal  processing  and  VLSI  array 
processing  (especiab^  SAR  processing). 

2.  the  design  methodology  of  VLSI  array  processors  for  these  applications. 

•  TNO-FEL  (applied  research  and  development)  and  the  Delft  University  of  Technol¬ 
ogy  (basic  research)  combine  their  expertises  on  the  fields  of  signal  processing,  signal 
processing  applications  (radar)  and  VLSI  array  processing.  The  leading  position  of 
TNO-FEL  in  the  field  of  SAR  processing  and  related  fields,  and  the  leading  position 
of  the  Delft  University  of  Technology  in  the  field  of  signal  processing  and  VLSI  array 
processing  are  maintained. 

•  The  research  is  presented  in  leading  signal  processing  and  radar  magazines  and 
conferences. 

•  A  PhD  thesis. 


3.3  Conditions  and  requirements 


Functional  requirements. 


Airborne:  The  SAR  processor  will  be  used  for  an  airborne  SAR. 

Real-time:  After  receiving  the  data  it  must  be  processed  and  displayed  instantaneously. 
On-board:  Processing  is  done  on  the  mrcraft. 

Operational  requirements. 

Image  quality:  The  a  SAR  processor  must  display  an  image  that  has  at  least  quick-look 
quality,  which  means  that  image  quality  and  resolution  are  of  less  importance.  A 
global  view  of  the  illuminated  area  is  displayed.  It  must  be  possible  to  check  the 
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equipment  on-line  whether  it  functions  well  or  not  and  to  judge  the  SAR  data  on 
its  usefulness.  The  image  quality  norm  may  be  increased,  as  long  as  the  data  is 
real-time  processed. 

Motion  compensation  (optional ):  The  quality  of  the  image  can  be  improved  by  phase  error 
compensation  with  motion  data  of  the  antenna  and  the  aircraft.  Motion  compensa¬ 
tion  is  implemented  if  real-time  processing  is  possible  and  if  quality  improvement  is 
necessary. 

Autofocussing  (optional):  The  quality  of  the  image  can  also  be  improved  by  autofocussing. 
The  conditions  for  implementation  are  the  same  as  for  motion  compensation. 

Design  requirements. 

VLSI  array  processing:  The  algorithms  for  SAR  processing  must  provide  a  high  through¬ 
put  of  SAR  data.  To  achieve  this,  the  algorithms  are  parallelized  as  much  as  possible 
and  mapped  on  array  processor  architectures.  An  ideal  real-time  on-board  SAR  pro¬ 
cessor  will  be  small,  which  implies  a  realization  in  VLSI  array  processors. 

Application  specific  IC’s  (ASIC’s):  To  meet  the  previous  requirement  a  choice  is  made 
for  an  application  specific  VLSI  implementation.  In  this  way,  speed  and  complexity 
can  be  optimized  as  much  as  possible. 


3.4  Project  outline 


The  project  consists  of  the  definition,  design  and  realization  of  a  real-time  SAR  processor. 
The  main  result  of  the  design  of  the  real-time  SAR  processor  is  an  operational  processor. 
The  first  stages  of  the  design  trajectory  shall  be  much  the  same  as  the  first  stages  of  the 
research  trajectory.  Basic  algorithm  and  architecture  design  is  required  to  come  to  this 
result.  The  architecture  d«ign  stage  results  in  a  theoretical  real-time  SAR  processor. 
In  the  first  stages  also  a  statement  about  the  feasibility  of  motion  compensation  and 
autofocussing  is  made,  based  on  research.  The  next  8t^;e8  wQl  lead  to  an  implementation 
and  a  realization  of  the  SAR  processor.  The  activities  are  focussed  on  the  application. 
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A  lay-out  is  designed  and  the  SAR  processor  is  manufactured.  The  processor  is  tested 
and  documented.  The  last  stages  differs  from  the  research  trajectory  since  theory  is  less 
important  in  the  realization  of  an  operational  product.  The  design  trajectory  of  the  real¬ 
time  SAR  processor  is  subdivided  in  the  following  stages: 

1.  Algorithm  design 

2.  Architecture  design 

3.  Lay-out  design 

4.  Implementation  and  realization 

5.  Testing 

6.  Documentation 

The  research  is  focussed  madnly  on  algorithm  design  and  architecture  design.  The  theo¬ 
retical  results  will  be  published.  The  intention  is  that  this  will  lead  to  an  PhD  graduation. 
A  detailed  description  of  the  project  stages  is  found  in  the  next  two  chapters.  In  chapter 
6  the  research  and  PhD  graduation  are  described  in  detail. 
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4  PROJECT  STAGES 

4.1  Algorithm  design 

Activities.  Algorithms  for  SAR  processing,  motion  compensation  and  autofocussing 
are  studied  extensively.  Experience  is  obtained  in  the  field  of  SAR  processing.  Useful 
algorithms  are  selected.  The  following  stages  are  distinguished: 

Literature  study:  Information  about  conventional  SAR  processing,  motion  compensa¬ 
tion  and  autofocussing  is  collected.  Several  digital  signal  processing  algorithms  are 
studied. 

Analysis:  The  echo  signal  including  phase  errors  is  analyzed.  A  matrix  description  of  the 
signal  is  derived.  The  echo  signal  is  of  a  convolution  of  the  transmitted  signal  (in¬ 
cluding  atmospherical  loss,  antenna  pattern,  etc.)  and  the  reflectivity  function  of  the 
illuminated  surface.  The  reflectivity  function  contains  the  geometrical  information 
of  the  surface.  * 

Synthesis:  An  inverse  operation  is  derived  to  restore  the  information  of  the  earths  surface. 
A  parallel  algorithm  is  determined  that  performs  this  algorithm  efficiently. 

Tools.  The  theoretical  background  to  perform  the  activities  mainly  consists  of  synthetic 
aperture  radar  theory  and  linear  algebra.  Software  tools  for  simulation  and  verification  of 
algorithms  are  MatLab,  C-programming  language  and  other  simulation  tools.  Hardware 
tools  are  a  SUN  Sparc2  station,  personal  computers  and  a  Convex  supercomputer.  The 
algorithms  are  tested  on  SAR  data  from  the  PHARS. 

Results.  The  results  of  this  stage  will  be  a  set  of  parallel  algorithms  to  perform  real-time 
SAR  processing,  which  can  be  mapped  on  parallel  VLSI  array  architectures. 
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4.2  Architecture  design 

Activities.  The  main  activity  is  the  mapping  of  the  algorithms  from  the  first  stage  onto 
array  architectures.  This  is  done  in  different  stages  [17]; 

Dependence  graph  (DG)  design:  The  algorithms  of  the  previous  stage  are  graphical  rep¬ 
resented.  The  graph  shows  the  dependence  of  the  computations  that  occur  in  the 
algorithms. 

Mapping  the  DG  onto  an  array  architecture:  The  regularity  of  the  DG  is  investigated 
and  the  best  mapping  methodology  is  determined.  One  or  more  array  processors  are 
designed. 

Tools.  The  behavior  of  the  architecture  design  is  simulated  on  different  levels  with  high 
level  simulation  languages.  It  is  tested  on  the  SAR  data  from  the  PHARS. 

Results.  The  architecture  design  stage  results  in  a  array  processor  architectures.  Their 
behavior  and  interconnections  are  described  in  detail.  The  architecture  is  specified  in 
processing  units,  control  units,  data  busses,  etc.. 


4.3 


Lay-out  design 


Activities.  The  processing  units  are  investigated  on  feasibility.  A  relation  between  ac¬ 
curacy  and  SAR  resolution  on  the  one  side  and  and  processing  speed  and  circuit  area  on 
the  other  side  is  determined.  An  optimum  is  determined  which  satisfies  the  requirements. 
The  arithmetic,  datacommnnication,  types  of  computations,  control  methodology,  etc.  are 
specified  from  the  requirements.  Before  the  actual  hardware  design  is  started  the  imple¬ 
mentation  of  proc^ng  units,  memory  units,  contrd  units,  etc.  are  specified.  Alternative 
implementations  are  for  example  ASIC’s  (Application  Specific  Integrated  Circuits),  trans¬ 
puters,  general  purpose  DSP’s.  After  this  specification  the  actual  hardware  ccHnpcments 
are  designed.  ASIC  dedgn  results  in  a  lay-out,  transput^  and  general  purpose  DSP  de^;n 
result  in  programming  code. 


Tools.  The  components  and  the  system  are  designed  with  hardware  simulation  languages. 
Lay-out  is  designed  with  (graphical)  lay-out  design  systems. 

Results.  The  result  is  a  set  of  hardware  components,  specified  and  designed,  and  a  system 
design.  The  complete  design  is  tested  on  SAR  data. 

4.4  Implementation  and  realization 

Some  hardware  components  like  general  purpose  DSP’s  and  transputers  are  bought  from 
external  suppliers.  The  programming  is  carried  out  at  TNO-FEL.  ASIC’s  are  manufac¬ 
tured  externally.  The  whole  system  is  assembled  at  TNO-FEL.  During  the  previous  stages 
this  stage  will  be  specified  in  detail. 


4.5  Testing 

Each  design  stage  has  its  own  test  and  verification  run.  The  designs  are  tested  off-line  on 
SAR  data.  This  raw  SAR  data  has  been  recorded  during  testflights  with  the  PHARS.  The 
several  components  of  the  system  are  tested  separately.  The  complete  system  is  tested  off¬ 
line  on  SAR  data.  At  the  moment  the  coUected  SAR  data  from  the  PHARS  is  recorded 
on-line  on  a  high-speed  taperecorder.  To  test  the  real-time  SAR  system  on  board  the 
SAR  data  will  also  be  sent  on-board  to  the  SAR  processor  for  immediate  processing  and 
displaying. 

4.6  Documentation 

Each  stage  is  documented  in  a  FEL-rep<^.  At  the  end  of  the  docnmntatkm  stage  a 
coherent  document  is  produced.  This  document  describes  the  design  trajectcffy,  the  system 
and  its  spedficatkms  and  the  user  speciilcatkms. 
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5  PROJECT  MANAGEMENT 

5.1  Time  planning 


The  time  planning  of  the  project  is  shown  in  figure  5.1.  Because  algorithm  and  architecture 
design  are  closely  related,  a  large  overlap  in  time  between  the  the  two  stages  is  reserved. 
Additional  to  the  project  stages  time  is  reserved  for  the  writing  of  the  PhD  thesis.  Parallel 
to  the  project  stages,  time  is  reserved  for  literature  study  and  writing  of  publications. 
These  activities  can  be  considered  as  part  of  the  subsequent  stages. 


Definition 
Literature  study 
Algorithm  design 
Architecture  design 

Lay-out  design 

Implementation  and 
Realization 

Testing 

Documentation 

Thesis 

Publications 


1995 


Figure  5.1:  Project  planning. 


5.2  Quality 

Two  qaafi^  norma  am  defined  to  meat  the  raquirementa.  The  first  deak  with  the  hardwam, 
the  aeomd  edth  the  imafa: 
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Hardware  quality.  In  order  to  perform  real-time  on-board  SAR  processing  the  hardware 
must  be  capable  of  processing  the  SAR  data  efficiently  and  fast.  This  results  in  the 
following  issues  during  the  design  stages: 

1.  The  algorithms  must  be  regular  and  parallel  and  have  low  computational  com¬ 
plexity 

2.  The  architectures  must  be  parallel,  pipelined  and  dedicated  to  array  processing 

3.  The  lay-out  must  be  dedicated  to  VLSI  array  processing 

Image  quality.  The  image  must  have  at  least  quick-look  quality,  which  means  that  image 
quality  and  resolution  are  of  less  importance.  A  global  view  of  the  illuminated  area 
is  displayed.  It  must  be  possible  to  check  the  equipment  on-line  whether  it  functions 
weU  or  not  and  to  judge  the  SAR  data  on  its  usefulness.  The  image  quality  norm 
may  be  increased,  as  long  as  the  data  is  real-time  processed. 


5.3  Communication  and  information 

The  scheme  in  figure  5.2  describes  the  information  flow  in  terms  of  input  and  output  per 
stage  (the  output  information  that  is  used  as  input  for  the  next  stage  is  marked  by  dotted 
arrows).  The  input  and  output  items  shall  be  described  briefly: 

Requirements:  The  requirements  are  defined  in  advance  (chapter  3). 

Theory.  The  theory  describes  why  a  particular  algorithm,  architecture  or  lay-out  is  cho¬ 
sen. 


Methodology.  The  methodology  describes  how  a  particalar  algorithm,  architecture  or 
lay-out  is  designed. 

Simuhted  behavior.  The  simulated  behavior  is  the  result  ol  the  nmulaUoa  runs. 

Spee^kationsr.  The  spedfleations  ate  derived  frena  the  simulated  bdmdor  aad  the  phjrs- 
kal  eoBstraiats,  i.e.  th«y  ate  updated  after  each  stage. 
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Bill  of  materials:  The  bill  of  materials  is  a  list  of  all  necessary  components  (including 
software  modules). 

Assembly  instructions:  The  assmebly  instructions  give  the  necessary  information  to  as¬ 
semble  the  processor. 

Expected  behavior.  From  the  specifications  of  the  software  and  hardware  the  expected 
behavior  of  the  processor  can  be  derived. 

Test  instructions:  The  test  instructions  give  the  necessary  information  to  test  the  pro¬ 
cessor. 
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Figue  5^:  lalormatioa  iowt. 


Page 

35 


Measured  behavior  and  specifications:  Erom  the  testruns  the  actual  behavior  and  specifi¬ 
cations  are  obtained. 

The  remaining  items  are  trivial.  After  each  stage  the  output  information  is  documented. 
In  the  documentation  stage  the  information  of  all  stages  is  collected  and  documented. 
Project  progress  is  reported  to  the  project  manager  (Ir.  H.  van  Es),  and  research  progress 
to  the  promotor  (Prof.Dr.Ir.  P.Dewilde). 


5.4  Organization 


The  project  is  carried  out  in  service  of  TNO  Physics  and  Electronics  Laboratory  (TNO- 
FEL).  The  research  is  carried  out  in  corporation  with  the  section  Network  Theory  of  the 
Delft  University  of  Technology  (TUD).  The  organization  of  the  project  is  as  foUows: 

Project  manager.  Ir.  H.  van  Es  (TNO-FEL) 

Promotor.  Prof.Dr.Ir.  P.  Dewilde  I(TUD) 

Mentor.  Dr.Ir.  E.  Deprettere  (TUD). 

Project  executor.  Ir.  L.  Bierens  (TNO-FEL). 
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6  RESEARCH 

6.1  Research  area 


As  mentioned  in  chapter  3  the  goal  of  this  project  is  to  come  to  a  deiinition  of  processor 
for  high  performance  real-time  on-board  SAR  processing.  Due  to  the  great  number  of  data 
and  the  high  computational  complexity  conventional  algorithms  are  inadequate  to  perform 
this  kind  of  SAR  processing.  The  mapping  of  these  algorithms  on  new  architectures  might 
increase  the  performance  but  not  sufficiently  enough  to  meet  the  demands.  Therefore  the 
first  part  of  the  research  of  this  project  is  focussed  on  algorithm  design.  The  second  part 
of  the  research  is  focussed  on  mapping  the  algorithms  on  new  architectures.  The  theory 
for  the  algorithm  design  is  found  in  the  fields  of  synthetic  aperture  radar,  linear  algebra 
and  digital  signal  processing.  Many  algorithms  for  digital  signal  processing  applications 
are  derived  with  linear  algebra  techniques  such  as  matrix  computations.  However,  this 
approach  is  rare  in  SAR  processing.  The  advantage  of  the  use  fineax  algebra  is  that 
algorithms  are  forced  to  be  highly  regular.  These  algorithms  can  be  mapped  on  high 
performance  array  architectures  via  a  systematical  design  methodologies  [17]. 


6.2  Research  description 


The  research  will  be  carried  out  mainly  in  the  algorithm  design  and  architecture  design 
stages.  In  the  last  year  of  the  project  additional  time  is  reserved  for  the  writing  of  the 
PhD  thesis. 

Algorithm  design.  In  the  algorithm  design  stage  the  following  research  activities  can 
be  distinguished: 

1.  As  preparation  for  the  actual  research  a  literature  study  to  matrix  computations, 
array  procesring,  SAR  processiiig,  etc.  is  carried  out. 

2.  Ctmventkmal  approaches  assume  that  transformation  of  the  SAR  signal  to  another 
d<nnain  (for  example  FFT  domain)  is  necessary  to  perform  operations  as  pulse  and 
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azimuth  compression.  With  linear  algebra  techniques  the  transformation  may  be 
avoided  and  processing  can  be  done  in  time  domain.  To  investigate  the  possibilities 
of  time  domain  processing  an  analysis  of  the  raw  SAR  signal  is  required.  I 

3.  After  the  raw  SAR  data  is  described,  SAR  processing  algorithms  are  designed.  The 
requirements  for  these  algorithms  are  regularity,  parallelism  and  low  computational 
complexity. 

Architecture  design.  In  the  architecture  design  stage  the  algorithms  from  the  previous 
stage  are  mapped  on  array  architectures.  This  requires  a  detailed  analysis  of  dependencies 
of  computations  that  occur  in  the  algorithm.  The  dependence  graphs  are  mapped  on  signal 
flow  graphs.  Finally  the  signal  flow  graph  are  mapped  on  array  processor. 

PhD  thesis.  The  intention  is  that  the  research  leads  to  a  PhD  graduation.  Therefore 
additional  time  is  reserved  for  the  writing  of  a  PhD  thesis. 

Parallel  to  the  activities  in  each  stage,  there  will  be  secondary  research  activities.  These 
activities  are: 


•  a  Uterature  study  to  keep  knowledge  up-to-date  (i.e  continuation  of  the  literature 
study  from  the  first  stage). 

•  the  writing  and  preparation  of  publications,  for  example  in  digital  signal  processing 
and  remote  sensing  magazines  and  on  conferences. 


6.3  Organization  of  the  PhD  graduation 


The  research  is  carried  ont  in  cooperation  with  the  section  Network  Theory  of  the  Ddft 
University  of  Tbchndogy  (TUD).  Prof.Dr.Ir.  P.  Dewilde  will  act  as  promotor  and  Dr.lr. 
E.Deprettere  wtQ  act  as  mmitor.  Formal  reporting  will  be  done  once  in  two  months  globally. 
Informal  discussions  will  take  place  more  frequently. 


4 


TNO  report 


BIBLIOGRAPHY 

[1]  D.A.  Ausherman,  A.  Kozma,  J.L.  Walker,  H.M.  Jones,  and  E.C.  Poggio.  Develop¬ 
ments  in  radar  imaging.  IEEE  Transactions  on  Aerospace  and  Electronic  Systems, 
20(4):363-400,  July  1984. 

[2]  D.  Blacknell,  A.  Freeman,  S.  Quegan,  LA.  Ward,  I.P.  Finley,  C.J.  Oliver,  R.G.  White, 
and  J.W.  Wood.  Geometric  accuracy  in  airborne  sar  images.  IEEE  Transactions  on 
Aerospace  and  Electronic  Systems,  AES-25(2);241-256,  March  1989. 

[3]  D.  Blacknell,  A.F.  Freeman,  R.G.  White,  and  J.W.  Wood.  The  prediction  of  geometric 
distortion  in  airborne  synthetic  aperture  radr  imagery  from  autofocus  measurements. 
IEEE  Transactions  on  Geoscience  and  Remote  Sensing,  GE-25(6);775-782,  November 
1987. 

[4]  D.  Blacknell,  LA.  Ward,  and  A.  Freeman.  Motion  compensation  and  geometric  distor¬ 
tion  in  airborne  sar  imagery.  In  Proceedings  ISPRS  Symposium,  Progress  in  Imaging 
Sensors,  pages  538-548,  September  1986. 

[5]  W.D.  Brown  and  D.C.  Ghiglia.  Some  methods  for  reducing  propagation-induced  phase 
errors  in  coherent  imaging  systems.  1.  formalism.  Journal  of  the  Optical  Society  of 
America  A,  5(6):924-941,  June  1988. 

[6]  A.  Di  Cenzo.  A  new  look  at  nonseperable  synthetic  aperture  radar  processing.  IEEE 
Transactions  on  Aerospace  and  Electronic  Systems,  24(3):218-224,  May  1988. 

[7]  J.  Dali.  Fast  convolution  algorithms  for  sar  processing.  In  Proceedings  International 
Conference  on  Radar  ’89,  pages  226-231,  April  1989. 

[8]  P.H.  Eichel,  D.C.  Ghiglia,  and  C.V.  Jakowattz  Jr.  Speckle  processing  method  for 
synthetic  aperture  radar  phase  correction.  Optics  Letters,  14(1):  1-3,  Januari  1989. 

[9]  R.  FVtbrizio.  A  high-speed  distal  processor  for  real-time  synthetic  aperture  radar 
imapng.  In  Proc.  IGARSS’87,  pages  1099-1103,  May  1987. 

[10]  LP.  Finley  and  J.W.  Wood.  An  investigation  of  synthetic  aperture  radar  antofocns. 
Internal  memorandnm,  rsre  memorandnm  3790,  Royal  Signals  k  Radar  EstaUish- 
mrat,  April  1989. 


;  I* 


TNO  report 


Page 

39 


[11]  G.  Franceschetti  and  G.  Schrinzi.  A  sai  processor  based  on  two-dimensional  fft  codes. 
IEEE  Transactions  on  Aerospace  and  Electronic  Systems,  26(2):356-366,  March  1990. 

[12]  D.C.  Ghiglia  and  D.C.  Brown.  Some  methods  for  reducing  propagation-induced  phase 
errors  in  coherent  imaging  systems.  2.  numerical  results.  Journal  of  the  Optical  Society 
of  America  A,  5(6):942-957,  June  1988. 

[13]  S.  Holm  and  A.  Maoy.  Synthetic  aperture  radar  processing  facility  based  on  a  parallel 
supercomputer.  In  Proc.  IEEE  International  Conference  on  Acoustics,  Speech  and 
Signal  Processing,  May  1989. 

[14]  P.  Hoogeboom  and  P.J.  Koomen.  The  pharus  project:  Design  of  a  polarimetric  c-band 
sar  in  the  netherlands.  In  Proceedings  AGARD,  May  1989. 

[15]  S.A.  Hovanessian.  Introduction  to  Synthetic  Array  and  Antenna  Imaging  Radars. 
Arthech  House,  1980. 

[16]  J.C.  Kirk.  A  discussion  of  Uigital  processing  in  synthetic  aperture  radar.  IEEE 
Transactions  on  Aerospace  and  Electronics  Systems,  AES-ll(3):326-337,  May  1975. 

[17]  S.Y.  Kung.  VLSI  Array  Processors.  Prentice  Hall,  1988. 

[18]  C.J.  Oliver.  Review  article:  Synthetic  aperture  radar  ims^ng.  Journal  of  Applied 
Physics,  22(4):871-890, 1989. 

[19]  M.P.G.  Otten.  Phase  error  estimation  for  pharus:  A  review  of  autofocus  and  related 
algorithms.  Internal  report,  no.  fel-89-c243,  TNO  Physics  and  Electronics  Laboratory, 
September  1989. 

[20]  M.P.G.  Otten.  Sax  processing  methods.  Internal  report,  no.  fel  1989-131,  TNO  Physics 
and  Electronics  Laboratory,  April  1989. 

[21]  M.P.G.  Otten.  Comparrison  of  sar  autofocus  algorithms.  In  Proceedings  Military 
Microroave  Conferenrx,  June  1990. 

[22}  C.  Pratio.  Autofocussing  synthetic  aperture  radar  images.  In  Proceedings  Interna¬ 
tional  Conference  on  Radar  ’89,  pages  314-319,  April  1989. 

[23]  K.W.  Przytula  and  J.G.  Nash.  Paralld  implementation  of  synthetic  aperture  radar 
alg(»ithms.  Journal  of  VLSI  Signal  Processing,  1:45-56, 1989. 


TNO  report 


Page 

40 


[24]  K.W.  Przytula,  J.G.  Nash,  and  S.  Hansen.  Fast  fourier  transform  algorithm  for  two- 
dimensional  array  of  processors.  In  Advanced  Algorithms  and  Architectures  for  Signal 
Processing  II,  pages  186-198,  1987.  Proc.  SPIE  826. 

[25]  S.  Quegan.  Measurement  of  geometric  distortion  in  airborne  sar  images.  In  Proceed¬ 
ings  IGARS’S4  Symposium,  ps^es  595-599, 1984. 

[26]  M.I.  Skolnik,  editor.  Radar  Handbook.  McGraw-Hill,  1970. 

[27]  A.M.  Smith.  A  new  approach  to  rangc-doppler  sar  processing.  International  Journal 
of  Remote  Sensing,  12(2):235-251, 1991. 

[28]  T.K.  Truong,  I.S.  Reed,  R.G.  Lipes,  A.L.  Rubin,  and  S.A.  Butman.  Digital  sar 
processing  using  a  fast  polynomial  transform.  IEEE  Transactions  on  Acoustics,  Speech 
and  Signal  Processing,  ASSP-22(2):419-425,  April  1984. 

[29]  F.T.  Ulaby,  R.K.  Moore,  and  A.K.  Fung.  Microwave  Remote  Sensing,  Active  and 
Passive,  volume  2.  Addison- Wesley,  1982. 

[30]  J.W.  Wood.  The  removal  of  azimuth  distortion  in  synthetic  aperture  images.  Inter¬ 
national  Journal  of  Remote  Sensing,  9(6):1097-1107, 1988. 

[31]  C.  Wu,  R.  Barkan,  W.J.  Karplus,  and  D.  Caswell.  Seasat  synthetic  aperture  radar 
data  reduction  using  parallel  programmable  array  processors.  IEEE  Transactions  on 
Geoscience  and  Remote  Sensing,  GE-20(3):352-358,  July  1982. 


Group  Leader 
R.IL  van  Es 


L.H.J.  Bierens 


REPORT  DOCUMENTATION  PAGE 


(MOD-NL) 


1.  DEFENSE  REPORT  NUKttER  (MOD-NL)  2.  RECIPIENTS  ACCESSION  NUMBER  3.  PERFORMING  ORGANIZATION  REPORT 

NUMBER 

TD91-3276  FEL-91-B295 


4.PRaiECT/rASK/WORKUNrrNO.  5.  CONTRACT  NUMBER 

22458 


6.  REPORT  DATE 
NOVEMBER  1991 


7.  NUMBER  OF  PAGES  8.  NUMBER  OF  REFERENCES  9.  TYPE  OF  REPORT  AND  DATES  COVERED 

40  (EXCL.  RDP  &  DISTR.UST)  31  RNAL 


10.  TITLE  AND  SUBTITLE 

A  REAL-TIME  SYNTHETIC  APERTURE  RADAR  PROCESSOR:  INTRODUCTION  AND  PROJECT  DESCRIPTION 


II.AUTHOR(S) 
L.H.J.  BIERENS 


12.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS<ES) 

TNO  PHYSICS  AND  ELECTRONICS  LABORATORY,  P.O.  BOX  96864, 2509  JG  THE  HAGUE 
OUDE  WAALSDORPERWEG  63,  THE  HAGUE,  THE  NETHERLANDS 


13.  SPONSORING4MONITORING  AGENCY  NAME(S) 

TNO  PHYSICS  AND  ELECTRONICS  LABORATORY,  THE  HAGUE,  THE  NETHERLANDS 


14.  SUPPLEMENTARY  NOTES 


15.  ABSTRACT  (MAXIMUM  200  WORDS,  1044  POSITIONS) 

THIS  REPORT  PRESENTS  A  PROJECTS  DESCRIPTION  FOR  THE  DESIGN  AND  REALIZATION  OF  A  REAL-TIME  SYNTHETIC 
APERTURE  RADAR  (SAR)  PROCESSOR.  THE  PRINCIPLES  OF  SAR  PROCESSING  AND  THE  STATE  OF  THE  ART  ARE 
DESCRIBED.  CONVENTIONAL  SAR  PROCESSING  TECHNIQUES  DEPEND  LARGELY  ON  FFT  TRANSFORMATION 
TECHNIQUES  AND  DEVICES.  THESE  OFFER  UMITED  POSSIBILITIES  TO  INCREASE  SPEED  AND  PERFORMANCE  OF  A  SAR 
PROCESSOR.  THEREFORE  THE  DESIGN  OF  NEW  ALGORITHMS  RATHER  THAN  ARCHITECTURES  OF  EXISTING 
ALGORITHMS  IS  SUGGESTED.  AFTER  THAT  NEW  ARCHITECTURES  WILL  BE  DESIGNED  DEDICATED  TO  THE  NEW 
ALGORITHMS.  BY  EXPLORING  PARALLEL  STRUCTURES  OF  THE  ALGORITHMS  AND  BY  USING  VLSI  ARRAY  PROCESSING, 
SPEED  AND  PERFORMANCE  CAN  BE  INCREASED  SIGNIFICANTLY.  THE  INTENTION  IS  THAT  THE  PROJECT  WILL  LEAD  TO 
A  PHD  GRADUATION. 


16.  DESCRIPTORS  IDENTIFIERS 

SYNTHETIC  APERTURE  RADAR  DIGITAL  SIGNAL  PROCESSORS 

SIGNAL  PROCESSING 

DESIGN 

ALGORITHMS 


17B.  BECMUTY  OABWICATION 
(OF  REPORT) 
UNOASSFIED 


ITBl  8ECURITY  CLASSmCATION 

(OFPA(M9 

UNClASaHEO 


17a  BECUmTYCLABBBTCATION 
(OFMnrnMCT) 

WCLASSnEO 


IB.  0«TRNRR10NMVAMEBJTY  BTAUMBfr 

UNLMTEO 


ITd.  tnURITY  OLASamCATION 
(OPirriM) 

UNCUWB) 


UNCLASSIPIB) 


Distribatielijst 


1.  Hoofddiiecteur  van  de  Hoofitgro^  Defeosieondeizodi:  TNO 

2.  Directeor  Wetenschappelijk  Oodeizodc  en  Qntwilckeling 

3.  HWO-KL 

4. 

t/m  HWO-KLu 

5. 

6.  HWO-KM 

7. 

t/m  Hoofd  TDCK 

9. 

10.  TUD,  ta.v.  prof.dr.ir.  P.  Dewilde 

11.  TUD,  t.a.v.  drJr.  E.  Depretter 

12.  Directie  FEL-TNO,  ta.v.  ir.  P.  Spohr 

13.  Direcde  FEL-TNO,  t.a.v.  dr.  J.W.  Maas,  daama  reserve 

14.  Archief  FEL-TNO,  in  bruikleen  aan  ir.  C.  Eberwijn 

15.  Archief  FEL-TNO,  in  bruikleen  aan  ir.  H.R.  van  Es 

16.  Archief  FEL-TNO,  in  bruikleen  aan  ir.M.P.G.Otten 

17.  Archief  FEL-TNO,  in  bruikleen  aan  ir.  P  J.  Kotnnen 

18.  Archief  FEL-TNO,  in  bruikleen  aan  ir.  P.  Hoogeboom 

19.  Archief  FEL-TNO,  in  bruikleen  aan  ir.  A.W.P.  van  Heijningen 

20.  Docummtatie  FEL-TNO 

21. 

t/m  Reserves 

25. 


iHdhwiwtilwKspiasMhatiaMIdaBiaslwrWaiwiMNwulp 


